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Abstract 

We propose to use the hadronic shape variables to test the assumption of "local 
duality" in the doubly-charmed B decays, B — * X C c S . We compute these shape 
variables at leading order in parton model, and discuss how they can be compared 
with the experimental values. 
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Recently there has been a lot of progress in understanding inclusive B decays us- 
ing the operator product expansion (OPE) and the heavy quark effective theory [1]. 
A better understanding of inclusive B decays will enable us to calculate the semilep- 
tonic and nonleptonic decay rates with less theoretical uncertainties, and therefore 
help us to get a more precise determination of the CKM matrix elements |V^| and 
\V u b\. In addition, the B semileptonic and nonleptonic decay rates also determine the 
semileptonic branching fraction of B decays 1 . 

The decay rates of B into inclusive semileptonic and nonleptonic final states have 
an OPE, of which the leading order reproduces the parton model result. However, 
the validity of the OPE is different for semileptonic and nonleptonic decays. For 
inclusive semileptonic B decays {b — > civ at quark level), the decay rate (and the 
differential spectra) can be expressed as an integral over some kinematic variable. In 
the complex plane of this kinematic variable, the only singularity is the physical cut on 
the real axis. It is possible to deform the integration contour for the decay rate (and 
the differential spectra) of this kinematic variable away from the physical cut in the 
complex plane. Along the deformed contour away from the physical cut, one hopes 
that the resonance effects associated with long distance hadronic production can be 
avoided, and perturbative QCD can be used reliably to compute the inclusive decay 
rate (and the differential spectra). This is called "global duality". It is generally 
believed to be valid. For inclusive nonleptonic B decays {b — > cud, b — > ccs at 
quark level), however, one has to invoke quark-hadron duality directly in the physical 
region in order to use the OPE (and perturbative QCD) to compute the decay rates. 
This assumption is called "local duality" . The validity of local duality relies on the 
energy release in the B decay being much larger than Aqcd, m which case there are 
many available states in the decay products. It is then likely that the decay rate is 
insensitive to the resonance effects associated with the production of particular hadron 
state, and one may apply local duality in the physical region to compute the B decay 

1 At one point it was thought that the theoretical value for the semileptonic branching fraction 
is too high to be consistent with the experimental value [2]. However, given the uncertainties 
associated with the theoretical prediction, it seems premature to talk about its discrepancy 
with the experimental value [3]. 
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rate. However, it is not clear whether the B meson is massive enough for local duality 
to apply. Especially in the case of the doubly-charmed B decays (6 — > ccs at quark 
level) , due to the large mass of the charm quark, the energy release A = — 2m c 
divided among three partons in the final state is not much larger than the QCD scale, 
and the resonance effects may be significant enough to destroy the validity of local 
duality. Therefore, it seems particularly important to test the validity of local duality 
for inclusive nonleptonic B decays in as many ways as possible, especially for the 
doubly-charmed decays. In this letter, we propose to use the shape variables of the 
doubly-charmed B decays as a test of local duality. 

The shape variables were first introduced in [4] for B decays into semileptonic 
final states and singly- charmed nonleptonic final states. In analogy we can define a 
similar shape variable for B decays into doubly-charmed final states {b — > ccs) 



T is the total decay width of B into doubly-charmed final states, m& is the b-quark 
mass. The first sum is over any hadronic pair for % ^ j, with and 
counted as different pairs. 9 is the angle between the two hadrons in pair. As 
noted in [4], the shape variable defined above is infrared safe, free from divergences 
due to collinear gluon and soft gluon emissions. Consequently one can use free quarks 
and gluons to calculate this shape variable, i.e., one can compute the differential decay 
rate at parton level, and sum over all possible parton pairs (and partons) to get d ^ s6 - 
Notice that for this observable there is no OPE, unlike for the inclusive decay rates. 
All that is being assumed is local duality. 

In the limit ^> m c ^> Aqcd, the c and the c in the decay products move along 
different directions with momenta much larger than Aqcd (i n b rest frame) in most 
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of the phase space. Because c(c) quark is heavy, the effect of long-distance soft QCD 
interaction on the direction of c(c) is suppressed. To leading order in 0(AQCD/ m c), 
the direction of c(c) quark is a physical observable. One can define a shape variable 
for production of charmed and anti-charmed final states only 

d£ c 1 v- /" ,r, \Pc\ Wc\ d 3 T 

% ' dE c dE^^———— , (2) 



d cos 9 T J mi, m b dE c dEcd cos 9 

where the summation is over all possible charmed and anti-charmed hadrons, 9 is the 
angle between the charmed hadron and the anti-charmed hadron. As noted above, 
one can invoke local duality to calculate this shape variable using the free quark 
picture, and it should be valid to leading order in 0(AQCD/ m c)- Notice that it is not 
clear at which order the QCD nonperturbative corrections arise for the parton model 
result of d ^ s0 [5]. Consequently the parton model result for d d ^e ought to be more 
reliable. 

The shape variables d d f s0 and ^osfl are computed in leading order in parton 
model for various values of m c /m]j. The results are shown in fig. 1. (We have taken 
m s = 0. d d ^ s g for m c = 0, which is for decay b uud, was first computed in [4].) 
Due to kinematics, the shape variables are all strongly peaked at cos 9 = — 1, and fall 
off rapidly as cos 9 goes from —1 to +1. The falloff is more steep for d d ^ , which 
implies that the hadrons containing c and c respectively in the decay products prefer 
to travel back to back. It will be interesting to see whether the experimental data 
reproduce this feature. The ^-function contribution to d ^ S Q at cos# = +1 is shown 
in Fig. 2 as a function of m c /mfo. 

Local duality has been used to compute these shape variables. To compare with 
the experimental measurements, one has to smear both the calculated and the mea- 
sured shape variables [6]. The particular details of the smearing function are not 
important, except for its width. Take d d ^e f° r example. Consider the contributions 
to it at hadronic level. There is a significant branching fraction for producing char- 
monium (cc) final states. Its contribution to xs a ^- function at cos# = +1. 
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This 5-function is not found in the free quark picture, as long-distance nonpertur- 
bative QCD effects are essential in forming cc bound states. There is also a large 
branching fraction for DD two-body final states, which give a 5-function contribu- 



tion at cos 6 = —1 in d d ^e • ^ n ^ ne ^ ree Q uar k picture d a ^o * s onr y strongly peaked at 
cos 6 — — 1. Therefore smearing is crucial in comparing the theoretically-calculated 
shapes with the experimentally-measured ones. 



It is useful to define the Fox- Wolfram moments [7] 
H(n) = 



(I a is — -/'„(<•(» t 



dcosO 



HJn) 



dcosO f c P n (cos( 
a cos 9 



(3) 



where P n (x) are Legendre polynomials. The first several moments are presented in 
Tables 1 and 2. Notice that H(l) = because of momentum conservation. In Table 1 
H(0) is subtracted from all the moments to remove the 5-function contribution to 
wfSg - Apart from the contribution from the 5-function in d ^ s g , the high moments 
go to zero asymptotically as n — > oo. These moments can be measured experimentally, 
which will test the applicability of local duality to B — > X C c S . 



n 


H(n) 


-H(0) 




me/mi, = 0.3 


rric/mi, = 0.37 


1 


-0.49 


-0.30 


2 


-0.22 


-0.13 


3 


-0.36 


-0.22 


4 


-0.29 


-0.17 


5 


-0.33 


-0.20 


6 


-0.30 


-0.18 



Table 1. Fox- Wolfram moments for -J^ts . 



4 



n 


if c (n) 


777 / 7 77 7 — X 


7T7 / 7T7 L — D ^7 


U 


5.8 x 10 


3.8 x 10 


1 
i 


-3.9 x 


-2.9 x 10" 2 


2 


2.4 x 1(T 2 


2.0 x 10~ 2 


3 


-1.4 x 1(T 2 


-1.3 x 10~ 2 


4 


8.7 x 10~ 3 


8.9 x 10~ 3 


5 


-5.5 x 1(T 3 


-6.0 x 10~ 3 



Table 2. Fox- Wolfram moments for d ^ 6 ■ 

New physics other than the Standard Model of electroweak interactions can have 
different Lorentz structures for the b decay amplitude, such as scalar-scalar type or 
(V + A) type. It turns out that these various forms of the decay amplitudes all 
produce the similar plots of the shape variables at leading order in parton model. For 
example, the first several Fox- Wolfram moments H c (n) from the (V + A) x (V — A) 
type Lorentz structure differ from the Standard Model results by less than 10% at 
leading order in parton model. So the shape variables are not sensitive probes for 
new physics. The experimental measurements of these shape variables will constitute 
a test of local duality, rather than a test for the possibility of new physics. 

In summary, we have computed the shape variables for the doubly-charmed inclu- 
sive B decays. The results found here for B decays also apply for B a and decays. 
Although we have only computed the leading order parton model contributions, it 
is straightforward to include the 0(a s ) perturbative contributions. Measurements of 
these shape variables will constitute a test of the applicability of local duality in the 
inclusive B decays. As there is no OPE for these shape variables, it is possible that 
local duality is more likely to fail for these shape variables than for the decay rates. 
It should be realized that understanding of the nonperturbative corrections to these 
shape variables is poor. If the experimental measurements of these shape variables 
corroborate the theoretical calculations, then it is possible that local duality holds in 
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the decay B — > A C c S despite the small energy release, and we will be more confident 
in using it to calculate the inclusive decay rate. 
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Fig. 1: Shape variables (a) d ^ s0 and (b) d ^e plotted over cos 9. The dashed line 
is for rric/mi) = 0, the solid line m c /m6 = 0.30, and the dot-dashed line rric/mi, = 0.37. 
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Fig. 2: The coefficient of 5(1 — cos 8) term in d ^ sS plotted over m c /m6. 
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